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We present a search for the decay B — > r v T in a sample of 88.9 x 10 6 BB pairs recorded 
with the BABAR detector at the SLAC _B-Factory. One of the two B mesons from the T(4S) is 
reconstructed in a hadronic or a semileptonic final state and the decay products of the other B in 
the event are analyzed for consistency with a B~ — > t~v t decay. We find no evidence of a signal and 
set an upper limit on the branching fraction of B(B~ — * t~v t ) < 4.2 x 1CP 4 at the 90% confidence 
level. 

PACS numbers: 13.20.He, 14.40.Nd, 14.60.Fg 



In the Standard Model (SM) the leptonic decay B~ — > 
T ~ v r Q proceeds via the annihilation of the b and u 
quarks into a virtual W boson. Its amplitude is thus 
proportional to the product of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element \V u b\ and the B me- 
son decay constant fg. The SM branching fraction is 
given by: 



"V T ) = 



G F m B 2 



1 



8tt 

(9.3 ±3.9) x 10 



rn 
nr 

-5 



f 2 B \Vub\ 2 TB 
(1) 



where Gf is the Fermi coupling constant, m T and tub 
are the r lepton and B~ meson masses, and tb is the 
B~ mean lifetime. We have used tb = (1.671±0.018) ps, 
\V ub \ = (3.67±0.47)xl0~ 4 , and Jb = (0.196±0.032) GeV 
(obtained from lattice QCD calculations) The 
branching fractions for e~V e and iiTv^ are helicity sup- 
pressed by factors of ~ 10~ 8 and ~ 10~ 3 , respectively. 
Physics beyond the SM, such as supersymmetry or two- 
Higgs doublet models, could enhance B(B~ — > r~z7 r ) 
by up to a factor of five through the introduction of a 
charged Higgs boson Q. 

A search for this decay is experimentally challenging 
due to the presence of at least two undetectable neutri- 
nos in the final state. No observation has been reported 
yet and the most stringent published limit on the decay 
is B(B~ -> t~i7 t ) < 5.7 x lfr 4 at the 90% confidence 
level |5J. 

The data used in this analysis were recorded with the 
BABAR detector at the PEP-II asymmetric e + e~ stor- 
age ring. The sample consists of 88.9 ±1.0 million BB 
pairs (81.9 fb" 1 ) collected at the T(45) resonance ("on- 
resonance ) and 9.6 fb 1 collected about 40MeV below 
the BB threshold ("off-resonance"). 

The BABAR detector is described in detail elsewhere . 
Detection of charged particles and measurement of their 
momenta are performed by a five-layer double-sided sil- 
icon vertex tracker and a 40-layer drift chamber, which 
operate in a 1.5-T solenoidal magnetic field. A detector 
of internally reflected Cherenkov light is used to identify 
charged kaons and pions. Photons and electrons are de- 
tected in an electromagnetic calorimeter consisting of an 
array of CsI(Tl) crystals. Muons and neutral hadrons are 
identified in the flux return, which is instrumented with 
multiple layers of resistive plate chambers. A GEANT4- 
based Q simulation of the BABAR detector, including 



machine backgrounds, is used to study signal event se- 
lection and background rejection. 

We first select a sample of events with one B-meson 
(the tag B) reconstructed in a hadronic or a semileptonic 
final state. The reconstruction constrains the kinematics 
and reduces the combinatorics in each event. This is 
critical since at least two neutrinos result from the B~ — > 
t~v t decay. All the neutral and charged particles not 
used for the tag B are assumed to come from the B- 
meson recoiling against it. We use two methods to search 
this recoil system for evidence of a B~ — > t~i7 t signal. 

In our first method, we reconstruct the tag B semilep- 
tonically. The semileptonic B-meson, B s \, is recon- 
structed as B + — > D°£ + iseX, where I = e,fi and X 
can be a 7, 7T°, or nothing. We select semileptonic 
B-decay events with several missing particles (such as 
neutrinos) by requiring at least one lepton with center- 
of-mass (CM) momentum (|p^|) above l.OGeV/c, zero 
event charge, a ratio of the Fox-Wolfram moments || 
H2/H0 < 0.9, and missing mass greater than 1.0 GeV/c 2 . 
Here, the missing mass is determined by subtracting the 
total energy and momentum of all reconstructed tracks 
and neutrals from the four-momentum of the T(45) sys- 
tem. We reconstruct D° mesons in the modes D° — > 
K + ir~, K + ir~ir~ir + , K + ir~7r°, and K®ir + ir~ and re- 
quire their reconstructed masses to be within three stan- 
dard deviations of the observed mean. The D° mesons 
are then paired with leptons with > l.OGeV/c to 
form D£ candidates. If the D° decay contains a charged 
kaon, the lepton must have the same charge as the kaon. 
The D° and lepton are required to originate from a com- 
mon vertex, but we do not mass-constrain the vertex fit. 
We assume that the only missing particle is a neutrino 
and calculate the cosine of the angle between the mo- 
mentum vectors of the D£ candidate and the B-meson, 



COSWr ne = , . Iz) 

The CM energy and momentum of the D£ candidate 
are E* Di and p^g, respectively. The B-meson energy is 
taken to be the beam energy, B^ oam . Calculated val- 
ues of cosOb.di may lie outside the physical range for 
events where the D£ candidate did not arise as pre- 
sumed, or due to detector energy and momentum resolu- 
tion. We place an asymmetric restriction on this variable, 
—2.5 < cos 6b,dz < 1-1) to admit D*° states where addi- 
tional decay products are present. If there is more than 
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efficiency correction from ali sources is determined to 
be 0.878 ± 0.076, and the corrected signal efficiency is 
(4.19±0.31 stat ±0.36 syst ) x lO" 4 . 

Probability density functions (p.d.f. 's) are constructed 
from the -Eextra distributions in signal (F (£" e xtra)s) and 
background (F(E ext ra)b) simulations. The £" C xtra dis- 
tribution for signal events is modeled as the sum of 
an exponential and two Gaussian distributions. The 
double- Gaussian models signal events where the X in 
B + — ► D°l + viX is a ir° or photon with a characteristic 
energy around 0.15 GeV. The exponential models signal 
events where such neutral particles are absent. To model 
background, as determined from Monte Carlo, we use a 
third-order polynomial. The p.d.f. 's are combined into 
an extended maximum likelihood function, 



FIG. 1: The distribution of -Eextra after applying all selec- 
tion criteria. The fit to the data and its components are also 
shown. The background is normalized to the data luminosity 
and the signal simulation is normalized arbitrarily. 



one acceptable Dl candidate, we choose the one whose 
D° mass is closest to the mean of the fitted distribution. 

After identifying the B s \, the remaining particles are 
required to be consistent with B~ — > t~v t where t~ — > 
e~v e v T or [i^v^Vt. Exactly one track with a small im- 
pact parameter relative to the primary vertex must re- 
main. The track must have p* < 1.2GeV/c, and must 
be identified as either an electron or muon. We reject 
e+e~ — * t+t~ events by restricting the angle of the track 
with respect to the event thrust axis (|cos^-^| < 0.9) 
and the minimum invariant mass constructable from any 
triplet of tracks in the event (M 3 min > 1.5GeV/c 2 ). 
In general, continuum events tend to peak sharply at 
| cos Op f\ — 1 and t + t~ events in particular tend to peak 
at values of M™ m below the r mass. 

The signal yield in the data is determined using the 
distribution of the total energy deposited in calorimeter 
clusters (with a minimum energy of 0.020 GeV) by neu- 
tral particles not associated with the D° decay in the 
semileptonic B s i candidate, -E oxtr a (Fig. P). This vari- 
able peaks near zero for signal while for background it 
rises with increasing -E ox tra- For -E e xtra < 1.0 GeV, we 
find from Monte Carlo simulations a signal efficiency of 
(4.77±0.35) x 10~ 4 and a background estimate of 124±7 
events. 

The signal efficiency quoted above is determined us- 
ing a detailed signal simulation. We study the differ- 
ences between simulation and data in the semileptonic B 
reconstruction, neutral-energy reconstruction, and lep- 
ton identification to derive an efficiency correction. The 
most significant effect comes from the B s \ reconstruction 
efficiency, and is determined using a sample of events 
in data and Monte Carlo simulations where both B 
mesons are reconstructed as B — > DlvX. The total 



C{s + b) 



j\\li,F{Ei) t + itbF[Ei) b 



(3) 



where Ei is the £" xtra i n the ith event, n is the total 
number of events in the data, and /j s and \ib are the signal 
and background yields to be fitted in the data. Studies of 
the choice of p.d.f. parameterization and of variations in 
shape suggest that the chosen p.d.f.'s yield a consistently 
conservative limit for the upper bound of the branching 
fraction. We fix the p.d.f. shape parameters and fit the 
data (Fig. The fit yields 14.8 ± 6.3 signal events and 
115.2 ± 11.8 background events. This signal yield has a 
statistical significance of 2.3c. 

We set a limit on the branching fraction at the 90% 
confidence level (C.L.) using the "CLs method" described 
in Refs. 00- We define our statistical estimator, Q, to 
be the fitted signal yield and compare the value of Q in 
data to its value in a large number of experiments gen- 
erated by sampling the likelihood function over a range 
of signal hypotheses. The uncertainty in the signal effi- 
ciency estimate is included by assuming a Gaussian un- 
certainty in the signal hypothesis. Using our fitted signal 
yield, efficiency, and the total number of B mesons in 
the data sample we determine that B(B~ — > t~i7 t ) < 
6.7 x 10~ 4 (90% C.L.). 

In our second method, we reconstruct the tag B can- 
didate, -Bhad, decaying into a set of purely hadronic final 

s*o 



The D* is reconstructed 
is a system of hadrons com- 



states, B+ _-> D^°X+. 
the mode D°ir°, and X^ 

posed of ni7r ± + n2K ± + 71371-° + n^K® where n\ = 1, ...5; 
ri2 = 0,1,2; 773 = 0,1,2; and 774 = 0,1. Rejection of 
background processes is based on two kinematic quan- 
tities: AE, the difference between the -Bhad and beam 
energies, and the beam-energy-substituted mass ttieSj 



m E s = \/[(s/2+p-p B ) 2 /E 2 } - \p B f 



(4) 



where ^/s is the total energy of the e + e~ system in 
the CM frame, and (E,p) and (Eb,pb) are the four- 
momenta of the e + e~ system and the -Bhad, respectively, 
both in the laboratory frame. 
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FIG. 2: Distribution of jties for the -Bhad candidates in data. 
The events lie in the region -0.1 < AE < 0.08 GeV. The 
solid curve shows the result of the fit with the sum of a Crystal 
Ball function (dashed curve) and an ARGUS function (dotted 
curve) . 



For each mode the toes distribution of the recon- 
structed candidates with —0.1 < AE < 0.08 GeV and 
toes > 5.21 GeV/ c 2 is fitted using the sum of a "Crys- 
tal Ball function" to model the signal component 
peaking at to_b and an "ARGUS function" to model 
the continuum and combinatorial B background. Fig- 
ure El shows the fit to the toes distribution for the 
Shad candidates in data. We define the signal region 
as -0.09 < AE < 0.06 GeV and ?7i ES > 5.27GeV/c 2 . 
We define a sideband region, 5.21 < toes < 5.26GeV/c 2 , 
to provide a control sample for studying continuum and 
combinatorial B background. The yield in the signal re- 
gion, as determined from the fit, is iVs had = (167.8 ± 
1.2 stat ± 3.0 sys t) x 10 3 . The error is dominated by sys- 
tematic uncertainty in the functional form of the peak at 
to b . 

We identify the r lepton using the following decay 
channels: r~ — *■ e~V e v Tl [i~v T T>^, n~v T , it~Tt°v T , and 
7r~7r + 7r~^v. We require the charged particles to be 
identified as leptons or pions, as appropriate. Mode- 
specific constraints are placed on the particles recoil- 
ing against the J5h a d- For the lepton and single-pion 
modes we reject events with 7r° or K® mesons in the 
recoil. The event is required to have zero charge and, 
in the recoil, at most one photon candidate not asso- 
ciated with a 7T°. Events with such a photon candi- 
date are accepted only if 50 < E 1 < 100 MeV (50 < 
Ej < 110 MeV for the r~ — > e~V e v Tl pTv T V^, and irv T 
modes) in the laboratory frame. Further requirements 
are made on the total missing momentum of the event, 
p m iss > 1.2GeV/c (> 1.4GeV/c for r~ -► vr-vr ^), the 
total momentum of the track(s) in the parent-S rest 



of two or three pions (0.60 < m W7r < 0.95 GeV/c 2 and 
1.10 < < 1.60GeV/c 2 for r" — > 7T 7T + 7T V T , 

0.50 < m„-„ < 1.00GeV/c 2 for t~ -► 7r - 7r i/ T )- 

We use detailed Monte Carlo simulations to deter- 
mine for each r decay channel the selection efficien- 
cies £i_weighted by the corresponding branching frac- 
tions 3j- The systematic uncertainties in selection ef- 
ficieny arise from tracking efficiency, neutral reconstruc- 
tion, particle identification, and ir° reconstruction. The 
total B~ — > t~i7 t selection efficiency (see Table UJ is 
(10.5 ± 0.2)%. Misreconstruction and contamination 
amongst the r-decay channels are taken into account. 

Continuum and combinatorial B background is deter- 
mined by extrapolating the ARGUS function from the 
toes sideband into the toes signal region. The back- 
ground that peaks in the toes signal region is determined 
from Monte Carlo simulations of B + B~ events. Events 
where a B is incorrectly reconstructed as a B + provide 
a negligible contribution. 

We correct the expected background, bi, to take into 
account possible dependencies of the fitted ARGUS 
shape on a given discriminating variable (p m i ss , invari- 
ant masses, etc.). The correction factor is the ratio of 
the background expectations determined using two sep- 
arate methods. In the first method, we estimate the 
background by scaling the number of events in the Toes 
sideband using the ARGUS signal-to-sideband ratio. In 
the second method, we bin each discriminating variable 
and reweight the number of events, bin-by-bin, using the 
ARGUS signal-to-sideband ratio for each bin. The sys- 
tematic error on bi is estimated as the deviation from 
unity of the total correction factor for each r-decay mode. 
The expected background and the total systematic un- 
certainty in each r-decay channel is reported in Table [i] 
along with the number rn of selected candidates in data. 

The systematic uncertainty in NB had (1.8%) is esti- 
mated as the change in the yield in the signal region in 
Fig. |21 when we use a double Gaussian as an alternative to 
the Crystal Ball function. Other models for the signal or 
the background distribution result in negligible changes. 

We observe a total of 15 B~ — ► t~v t candidates, 
which is consistent with the expected background of 
17.2 ± 2.1 sta t ± l-3 sys t events. The distribution of these 
events is also consistent with background. 

Wc determine the B~ — > r~ i7 T branching fraction from 
the number of signal candidates Si expected for each r 
decay mode, where = NB had B(B~ — ► T~i7 T )ei. The 
results for each decay channel are combined using the 
estimator, Q. Here we define Q to be C(s + b)/C(b), 
where 



»=i 



n , ! 



frame (p n - > 1.2GeV/c for r 
1.6GeV/c for r - — > 7r~7r + 7r~^ T ) 



> 



and the invariant mass 



(5) 

are the likelihood functions for signal-plus-background 
and background-only hypotheses and n c h is the total 
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TABLE I: Branching fraction (6)0! efficiency (e%), expected 
background (bi) with statistical and systematic errors, and 
observed data candidates (ru) for each reconstructed r decay 
mode. 



selection £>(%) Ei(%) bi rii 

~evv 17.84 ±0.06 3.4 ±0.1 0.7 ±0.4 ±0.1 2 

Hvv 17.37 ±0.06 1.9 ±0.1 0.9 ±0.5 ±0.1 

■ku 11.06 ±0.11 2.6 ±0.1 1.3 ±0.6 ±0.2 2 

■k~-k + -k~v 9.52 ±0.10 0.6 ±0.1 4.3 ±1.0 ±0.3 4 

7r~7T°z/ 25.41 ±0.14 2.0 ± 0.1 10.0 ± 1.6 ± 1.3 7 

all 81.20 ±0.22 10.5 ± 0.2 17.2 ±2.1 ±1.3 15 
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number of reconstructed r-decay channels. 

Since we have no evidence of signal we set an upper 
limit on the branching fraction. The statistical and sys- 
tematic uncertainties in the expected background are in- 
cluded in the estimator Q by convolving the likelihood 
functions with a Gaussian distribution having as stan- 
dard deviation the combined statistical and systematic 
errors in the background estimate 0|. We determine 
that B(B- -> t-i7 t ) < 4.2 x 10" 4 (90% C.L.). 

To combine the results from the statistically indepen- 
dent hadronic and semilcptonic samples, we first calcu- 
late the likelihood ratio estimator, Q = C(s + b)/£(b), us- 
ing the likelihood functions from each method. We create 
a combined estimator from the product of the semilcp- 
tonic (Qsi) and hadronic (Qhad) likelihood ratio estima- 
tors, Q = Qsi x Qhad- The measured branching fraction, 
which is the value that maximizes the likelihood ratio 
estimator, is (2.3^^3) x 10~ 4 . The lower one-standard- 
deviation bound does not include zero because of the 
small excess of signal events observed in the semilep- 
tonic analysis. Since this value is compatible with a zero 
branching fraction, we set a combined upper limit, 

B{B- -» t-v t ) < 4.2 x 1(T 4 (90% C.L.). (6) 

The semileptonic analysis does not contribute signifi- 
cantly to the combined limit because of the observed 
small excess of signal events. 

We use Eq.^ Eq.El and the measured value of \ V u b\ to 
set a limit on f B . We find f B < 0.510 GeV (90% C.L.). 

In conclusion, we have searched for B~~ — ► t~i7 t in 
the recoil of hadronic and semileptonic B decays. We 
have set the most stringent upper limit to date on this 
process. 
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